Thermal quadrupole method with internal heat sources

A new method based on the thermal quadrupoles technique for heat transfer modelling in multilayered slabs with heat sources is proposed. Classical thermal quadrupoles use hyperbolic functions and numerical problems occur according to the argument value that depends on thermophysical and geometrical properties as well as characteristics times. We propose a new formulation based on exponential function with negative argument. Using this formulation in the classical equivalent impedance network allows to compute efficiently the thermal behaviour of multilayered slabs with internal heat sources whatever the time and the thermophysical properties. This approach is applied in order to simulate heat transfer in three different multilayered materials with heat sources. These simulations show the capability of such a methodology to simulate time and space multiscale heat diffusion problems.